function [Hv] = func_H( a,b )
%Hy Возвращает вектор фокусоровки для нашей антенной системы
%   Входные данные - азимут  и угол места
[kx,ky,kz]=sph2cart(b,a,1);

k=[kx;ky;kz];

r{1}=[1/4;-1/4;0];
r{2}=[-1/4;-1/4;0];
r{3}=[-1/4;1/4;0];
r{4}=[1/4;1/4;0];

NAnt=length(r);


phi=nan(NAnt,1);

for i=1:NAnt
    
phi(i)=2*pi*(r{i})'*k;

end

Hv=exp(1i*phi);


end













